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The  integral-equation  method  developped  in  Part  I  is  applied  to  a  Pierce- 
type  diode  (Pierce  1944)  whose  external  circuit  involves  a  resistor,  an  in¬ 
ductor,  and  a  signal  generator.  The  general  linear  pertnrbational  problem 
is  solved  analytically  for  the  small-amplitnde  quantities  je{t)  (external- 
circuit  current  density)  and  &(x,  t)  (electrostatic  field).  Each  of  these 
quantities  can  be  constructed  from  a  "spatial*  Green’s  function  (describ¬ 
ing  initial  perturbations  of  the  plasma),  a  "temporal*  Green’s  function 
(describing  external-generator  signals),  and  two  functions  associated  with 
the  initial  state  of  the  external  circuit.  The  solutions  generally  exhibit  an 
initial  transient  and  an  asymptotic  part,  the  latter  being  a  superposition 
of  eigenmodes  only.  Systematic  numerical  results  for  eigenfrequencies  and 
eigenmode  profiles  in  some  typical  parameter  regions  demonstrate  that  the 
linear  response  and  stability  behavior  of  the  diode  system  may  substan¬ 
tially  depend  on  the  properties  of  the  external  circuit. 
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[.  INTRODUCTION  AND  SUMMARY 

-^Thia  paper  constitute*  Part  II  of  a  study  proposing  an  integral-equation 

approach  to  linear  longitudinal  oscillations  in  collisionless  plasma  diodes  with 

- - 

thin  sheaths.  The  method  was  formally  developed  in  Part  I*  and  is  comprehen¬ 
sive  in  that  it  (i)  handles  the  combined  initial-value  and  external-perturbation 
problem,  thus  including  the  eigenmode  problem  as  a  special  case,  and  (ii)  allows 
for  explicit  and  simultaneous  inclusion  of  three  basic  constituents  of  bounded 
plasma  systems,  namely  the  plasma  itself,  the  boundaries,  and  the  external  cir¬ 
cuits  As  discussed  in  Subsec.  m.  A/I, 3  The  present  restriction  to  one-dimensional 
diode  geometry  is  not  inherent  to  the  integral-equation  approach  as  such  and 

is  likely  to  become  dispensable  at  a  later  stage  of  development  of  the  theory. 

-f-hi  f 

In  the  present  paper,  the  method  developed  in  Part  I  m  employed  to  solve 
an  extended  Pierce-type  problem  involving  a  non-trivial  external  circuit.  Our 
goal  M^i)^ to  demonstrate  the  applicability  and  handling  of  the  method,  and 
show  in  numerical  detail  that  the  external  circuit,  an  essential  part  of 
any  bounded  plasma  system,  may  crucially  influence  the  overall  dynamic  and 
stability  behavior.  The  paper  proceeds  as  follows. 

Section  II  gives  a  description  of  the  Pierce  problem3-90  and  some  of  its 
implications.  FVom  the  discussion  and  literature  survey  in  Subsec.  II. A  it  is 
concluded  that  Pierce-type  instabilities  are  of  great  practical  importance,  e.g., 
in  various  laboratory  plasma  devices,  and  have  recently  attracted  considerable 
interest  in  the  context  of  such  modern  concepts  and  applications  as  double 
layers  or  inertial  confinement  fusion.  Areas  are  pointed  out  in  which  more  work 
is  needed,  among  them  the  complexes  of  external-circuit  effects  and  plasma- 
wall  interaction.  In  Subsec.  II. B  we  present  an  argument  according  to  which 
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the  Pierce  diode  represents,  under  certain  circumstances,  a  rough  but  still  useful 
first  approximation  to  the  one-emitter  plasma  diode  or  single-ended  Q  machine. 

In  Sec.  Ill,  the  method  developed  in  Part  I  is  formally  applied  to  derive 
the  integral  equations  governing  the  extended  Pierce-type  problem  considered, 

which  differs  from  the  “classical*  (Le.,  short-circuited)  one3  in  that  it  includes 

n 

an  external  circuit  j^ohring  a  signal  generator,  a  resistor,  and  an  inductor  as 
shown  in  Fig.  1.  Subsection  IRA  gives  the  formal  specifications  of  this  pertur- 
bational  problem,  i.e.,  the  level-one  quantities  defined  in  Table  I/I.  In  Subsec. 
m.B  we  try  to  foresee  the  structure  of  the  result  and  conclude  that  the  general 
solution  can  be  obtained  by  linear  superposition  of  four  different  types  of  “ba¬ 
sk  modes*,  two  of  whkh  are  identical  with  the  “spatial*  and  the  “temporal* 
Green’s  functions,  respectively.  The  level-two  through  level-four  quantities  of 
Table  I/I  are  calculated  in  Subsec.  IRC  and  used  in  Subsec.  HI.D  to  explkitly 
construct  the  integral  equations  (37) /I— (40) /I  as  well  as  the  equivalent  matrix 
equation  (48) /L  For  the  particular  case  under  consideration,  this  system  turns 
out  to  reduce  to  just  two  integral  equations,  the  unknowns  being  the  Laplace- 
transformed  perturbations  j«(u>)  (external-circuit  current  density)  and  &(z,u) 
(electroetatk  field). 

In  Subeec.  IV.  A,  these  integral  equations  are  solved  for  the  Laplace-  trans¬ 
formed  bask  modes,  whkh  are  Laplace  inverted  in  Subsec.  IV. B  into  the  time- 
dependent  bask  modes.  This  completes  the  general  solution  of  the  combined 
external- perturbation  and  initial-value  problem.  Each  of  the  basic  modes  ex¬ 
hibits  an  initial  transknt  as  well  as  an  asymptotic  part  and  can  be  represented 
as  a  linear  superposition  of  eigenmodes.  The  related  eigenffequencies  are  the 


seroe  of  the  characteristic  equation  (or  'dispersion  relation*),  whose  unnormal¬ 
ised  form  is  given  by  Eq.  (45). 

The  nondimensional  eigenmode  problem  is  dealt  with  analytically  in  Sec. 
V,  and  numerically  in  Sec.  VI.  In  Subsec.  V.A,  the  characteristic  equation  is 
re-written  in  nondimensional  form  and  found  to  depend  on  the  three  dimen¬ 
sionless  parameters  a  (“Pierce  parameter”),  p  (external-circuit  resistance),  and 
A  (external-circuit  inductance).  In  Subsecs.  V.B  and  V.C  we  discuss  some  an¬ 
alytic  properties  of  the  eigenfrequencies  and  eigenmode  profiles,  respectively. 

Throughout  Sec.  VI,  the  classical,  short-circuited  Pierce  diode3  (a  >  0,  p  = 
0,  A  *  0)  is  used  as  a  standard  reference  to  which  all  other  cases  are  conveniently 
compared.  It  is  reviewed  in  Subsec.  VI. A,  where  numerical  data  describing 
its  eigenfrequencies  and  eigenmode  structures  in  three  typical  a-intervals  are 
presented,  and  a  classification  scheme  for  the  eigenfrequencies  is  introduced. 
Subsection  VLB  deals  with  the  “purely  resistive*  case  (o  >  0,p  >  0,A  » 
0)  and  demonstrates  that  finite  external-circuit  resistance  can,  among  other 
things,  stabilise  a  classically  unstable  mode  and  change  the  character  of  the 
dominant  mode  from  oscillatory  to  non-oocillatory.  In  the  “purely  inductive* 
case  (a  >  0,p  =  0,  A  >  0),  which  is  considered  in  Subsec.  VI.C,  we  find  a  new 
oscillatory  mode  which  may  completely  change  the  linear  dynamic  behavior  and 
even  become  unstable  for  a- values  at  which  the  classical  Pierce  diode  is  stable. 
Some  general  aspects  and  implications  are  discussed  in  Subsec.  VI. D. 

Finally,  our  main  conclusions  are  summarised  in  Sec.  VII. 
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n.  THE  PIERCE  PROBLEM 


A.  Pierce  diode  and  Pierce  instability 

The  plasma  diode  model  given  by  Pierce  about  fourty  years  ago,3  which 
has  ever  since  been  reconsidered  in  various  modifications,4-30  is  a  particularly 
simple  special  case  of  the  collisionless  one-dimensional  configurations  described 
in  Subeec.  II.  A.  The  ions,  which  are  assumed  to  be  immobile,  provide  a  uniform 
neutralising  background  of  density  ftp,  immersed  between  the  left-hand  (x  =  0) 
and  right-hand  (x  =  L)  electrodes.  The  electrons  are  constantly  injected  at 
the  left-hand  electrode  as  a  cold  beam  with  density  ftp  and  velocity  0,  and  are 
absorbed  when  hitting  either  electrode.  Both  electrodes  are  invariably  kept  at 
the  same  potential,  which  corresponds  to  assuming  an  external  short-circuit. 

In  what  follows,  the  system  just  described  will  be  called  the  "classical" 
Pierce  diode.  Out  of  the  various  equilibrium  states  the  classical  Pierce  diode 
can  exhibit,  the  original  Pierce  paper3  was  only  concerned  with  the  one  in  which 
the  electrons  are  uniform  as  well.  This  special  case  will  be  referred  to  as  the 
"uniform  classical”  one. 

Nondimenskmal  analysis  shows  that  the  equilibrium  (d/dt  =  0)  and  dy¬ 
namic  states  of  the  classical  Pierce  diode  can  be  conveniently  classified  in  terms 
of  the  single  parameter 


where  wp  =  (4*r ftpe3/m)l  is  the  electron  plasma  frequency,  e  is  the  positive 
elementary  charge,  and  m  is  the  electron  mass.  The  parameter  a,  henceforth 
referred  to  as  "Pierce  parameter” ,  is  proportional  to  the  ratio  of  electron  transit 


time  to  electron  plasma  period  and  thus  provides  a  rough  measure  of  the  relative 
importance  of  plasma  versus  boundedness  effects. 

In  contrast  to  the  classical  Pierce  diode,  any  similar  plasma  configura¬ 
tion  characterised  by  a  strong  electron  current  between  two  electrodes  will  be 
referred  to  u  “Pierce-type’  or  “modified  Pierce*  diode.  Pierce-type  config¬ 
urations,  although  not  always  recognised  as  such,  are  widely  encountered  in 
practice,  cf.  below. 

Despite  its  simplicity  and  degeneracy,  the  uniform  classical  Pierce  diode 
is  of  fundamental  importance  because  it  qualitatively  retrieves  some  typical 
bounded-system  behavior  yet  is  still  tractable — to  some  extent,  at  least — with 
reasonable  mathematical  effort.  Most  importantly,  it  exhibits  the  “Pierce  insta¬ 
bility”,  which  destroys  the  electron  equilibrium  more  or  less  violently,  depending 
on  the  value  of  a.  Numerical  data  for  the  classical  Pierce  instability  may  be 
found,  e.g.,  in  Ref.  33,  Ref.  34,  and  Subsec.  VI.A  below. 

Since  the  classical  Pierce  instability  disappears  in  the  infinite-plasma  limit 
(Le.,  for  L  — ►  oo,  and  hence  a  — *  oo),  it  is  inherently  due  to  the  finite,  bounded 
character  of  the  system.  One  reason  for  this  striking  difference  between  fi¬ 
nite/bounded  and  infinite/unbounded  systems  can  be  spotted  in  the  different 
energetics:  For  a  shorter  diode,  the  energy  exchanged  with  the  electrodes  via 
particle  injection  and  absorption  is  of  greater  relative  importance,  so  that  the 
system  may  rapidly  evolve  to  a  final  state  whose  time-averaged  energy  is  signif¬ 
icantly  different  as  compared  to  the  initial  one;  if  the  diode  becomes  sufficiently 
long,  this  energy  exchange  is  relatively  unimportant  and  so  the  initial  energy  is 
practically  conserved.  A  second  important  difference  comes  about  through  the 
external  circuit,  whose  main  effect  is  to  control  the  surface  charges  sitting  on 


the  electrodes.  These  surface  charges  give  rise  to  a  vacuum  field  that  is  super¬ 
imposed  on  the  space-charge  field  of  the  plasma  and  is  inversely  proportional 
to  the  diode  length  L,  thus  becoming  less  effective  for  a  longer  diode. 

The  Pierce  instability  may,  generally  speaking,  occur  in  any  Pierce-type 
system  but  will  usually  not  be  as  violent  as  in  the  idealised  classical  case.  When¬ 
ever  it  appears  it  is  likely  to  play  a  non-negligible  or  even  dominant  role  in  the 
linear  and/or  nonlinear  dynamics  of  the  system.  It  may,  e.g.,  cause  a  rapid 
overall  transition  of  the  electron  and  potential  distributions  to  a  new,  “tempo¬ 
rary*  equilibrium  state,  which,  after  reaction  of  the  ions  on  their  much  longer 
time  scale,  can  result  in  double  layers, 2a.34.aT.3s.3s.42  phase-space  holes,30,36  or 
relaxation  oscillations.7'0,11,15'17'18,26,20,41,42  The  Pierce  instability  may  alter¬ 
natively  give  rise  to  electron  turbulence30,34,43  which,  in  turn,  can  influence  the 
ion  dynamics  via  the  ponderomotive  force.  In  all  cases,  the  Pierce  instability 
sets  an  upper  limit  to  the  time-averaged  electric  current  that  can  be  passed 
through  the  system. 

In  short,  the  Pierce  instability  is  not  just  of  academic  interest  but  rather 
requires  serious  consideration  whenever  one  is  dealing  with  bounded,  current- 
carrying,  weakly  collisional  plasma  systems.  It  has,  accordingly,  been  invoked  in 
the  context  of  various  practical  applications  and  devices  such  as  electron-beam 
generation,10,20  microwave  generation10,33,34  and  amplification,10  thermionic 
converters,7-0,11,16,18,26,20  drift  tubes,13  triode  lasers,10  or  ion-beam  neutral¬ 
isation  for  inertial-confinement  fusion.31,32,37,30  Its  relevance  to  low-density 
single-ended  Q  machines41,42  is  discussed  in  Subsec.  II.B  below.  Experimental 
evidence  of  the  Pierce  instability  can  be  found,  e.g.,  in  Refs.  5,  11,  22,  27,  and 
40-42. 
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Since  the  appearance  of  the  classical  Pierce  paper,  numerous  extensions 
and  modifications  of  the  Pierce  instability  have  been  studied  analytic- 

3^-9.8.10.12-14.16-23.28.27.29-37.39  ^  numerically  .7»1116-24.25.28.30-32.38 

These  studies  were  based  on  fluid,  3-6.10.12-14.1e-20.22.23.27.30-36.37.30 
kinetic, l*-*-31-36-3*-3*-38-3*-3*  or  particle7-8  11,16,34-30-33  models.  Beams  of  fi¬ 
nite  thermal  spread  were  considered  in  Refs.  6-9,  11,  15,  21,  23-26,  28-30, 
32,  36,  38,  and  39,  and  effects  of  finite  ion  mass  were  dealt  with  in  Refs.  4, 
7-13,  15,  17,  22-24,  27,  28,  31,  36,  and  38.  Apart  from  a  few  remarks  on 
collisional  effects  in  Ref.  10,  all  papers  mentioned  considered  the  collisionless 
case  only.  Most  studies  were  concerned  with  one-dimensional  geometry;  two- 
dimensional  effects  were  considered  in  Refs.  3,  5,  13,  20,  31,  and  32,  but  no 
three-dimensional  treatment  is  known  to  us.  Non-uniform  equilibria  were  ad¬ 
mitted  in  Refs.  6-9,  11,  12,  15,  18,  19,  21,  24-26,  28-31,  34,  and  38.  Both 
linear3-8-*' 10- 13’13,  i7-23.2e,27.29-33.3«, 37,39  and 

nonlinear7’8'11’14-1*-34'36-3*-30-33-34-3*  properties  of  the  Pierce  instability  have 
been  studied.  Except  for  some  considerations  on  electromagnetic  effects  in  Ref. 
13,  all  references  cited  employ  the  electrostatic  approximation. 

Given  the  fair  amount  of  literature  that  already  exists  on  the  Pierce  in¬ 
stability,  the  question  arises  as  to  what  remains  to  be  done.  In  our  opinion, 
there  are  two  basic  lines  along  which  progress  can — and  should — be  made,  (i) 
Several  areas  can  be  located  in  which,  to  our  knowledge,  very  little  detailed 
work  has  been  done.  E.g.,  the  role  of  non-trivial  external-circuit  elements, 
though  pointed  out  to  be  potentially  important,13  has  been  given  consideration 
by  only  very  few  authors.*'1318'37,3*  Another  example  is  the  fact  that  there 
does  not  seem  to  exist  any  treatment  including  realistic  plasma-wall  interaction 


processes  such  as  inelastic  reflection  or  secondary-particle  emission.  As  a  third 
example  we  wish  to  mention  that  only  little  analytic  work  has  been  devoted  to 
the  nonlinear  properties  of  the  Pierce  instability,141®-34-36  which,  however,  in¬ 
dicates  that  this  is  a  promising  area  where  much  more  progress  could  be  made, 
(ii)  In  view  of  the  complexity  of  realistic  systems,  it  seems  desirable  to  proceed 
towards  treatments  including  more  of  the  relevant  physics  simultaneously.  E.g., 
all  linear  analyses  available  are  concerned  with  the  eigenmode  problem — i.e., 
with  a  time  dependence  exp(at) — rather  than  with  the  full  initial-value  and 
external-excitation  problem. 

Generally  speaking,  we  believe  that  the  fundamental  importance  of  the 
Pierce  problem  requires  and  justifies  further  theoretical  efforts  towards  a  more 
complete  understanding.  Thus,  the  present  paper  contributes  to  two  aspects 
which,  according  to  the  above  literature  survey,  have  not  yet  been  treated  to 
a  significant  extent:  (i)  The  complete  initial-value  and  external-perturbation 
problem  is  solved  in  Secs.  IV  and  V,  and  (ii)  a  fair  amount  of  systematic  nu¬ 
merical  data  describing  the  influence  of  external-circuit  resistance  and  external- 
circuit  inductance  on  the  linear  eigenfrequencies  and  eigenmode  profiles  is  pre¬ 
sented  in  Sec.  VI.  The  results  of  particle  simulations  extending  into  the  nonlin¬ 
ear  regimes  of  the  Pierce  instability  will  be  presented  elsewhere.43 


B.  Cold-beam,  thin-sheath  approximation  of  the  positively  biased  one- 
emitter  plasma  diode  or  single-ended  Q  machine 

In  this  subsection  we  are  concerned  with  the  well-known  one-dimensional 
model  of  the  collisionless  one-emitter  plasma  diode  considered,  e.g.,  in  the  con¬ 
text  of  thermionic  converters44,46  or  low-density  single-ended  Q  machines.46-48 


Oar  goal  here  is  to  demonstrate  that  the  positive-bias  equilibrium  states  of  this 
model  can,  to  a  certain  extent,  be  approximated  by  the  uniform  classical  Pierce 
configuration.3  This  argument  is  of  practical  importance  because  it  allows  one 
to  obtain,  with  relative  ease,  quantitative  estimates  of  stability  domains  and 
linear  growth  rates.  Moreover,  the  numerical  results  following  from  this  model 
seem  to  indicate  that,  under  normal  operation  conditions,  single-ended  Q  mar 
chines  should  always  be  Pierce  unstable,  a  fact  which  apparently  has  not  been 
pointed  out  previously  but  could  possibly  be  of  relevance  in  explaining  the 
various  dynamic  phenomena  observed  experimentally.41-43  For  details  of  the 
model,  the  reader  is  referred  to  Refs.  44-48  and  many  others  cited  there.40 

The  equilibrium  states  in  question  are  characterised  by  the  potential  dis¬ 
tribution  types  1  (mono  tonic  ally  increasing)  and  2  (one-minimum)  as  defined  in 
Fig.  3  of  Ref.  46  and  Fig.  la  of  Ref.  4$.40  For  both  of  these  distribution  types, 
the  plasma  potential  Vv  is  positive  with  respect  to  that  of  the  emitter  surface, 
P(x  =  0)  =0.  Hence,  the  electrons — which  are  assumed  to  leave  the  emitter 
surface  with  a  half  Maxwellian  velocity  distribution — always  enter  the  plasma 
region  as  a  distinct  beam  with  a  truncated  Maxwellian  distribution.  The  cut¬ 
off  velocity,  average  velocity,  and  effective  temperature  of  this  distribution  are 
respectively  given  by4* 


r~i 

1 

<N 

II 

O 

(2o) 

/2  kT\*  exp  (— me3/2/cT) 

\  erfc  [oc(m/2*r)i] 

(26) 
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where  vm  is  the  minimum  potential  (cf.  Fig.  1/1;  for  a  monotonic  ally  increasing 
potential  distribution  as  shown  in  Fig.  1  we  have  Vm  =  0),  k  is  Boltsmann’s 
constant,  T  is  the  emitter  temperature,  and  erfc  is  the  complementary  error 
function.  For  a  ‘long*  system  such  as  the  Q  machine,  the  nondimensional  quan¬ 
tities  fe/f  =  Teff/T,  =  0e/v,  (where  v.  =  (2*r/m)i),  and  <fiav  =  0av/v,,  as 
well  as  several  others,  depend  only  on  the  “neutralisation  parameter*  &,  which 
is  defined  as  the  ratio  of  desorbing-ion  over  emitted-electron  density  at  the 
emitter  surface.4®'49  Figure  2b  of  Ref.  48  shows  that  at  &  -  0.405 — i.e.,  at  the 
transition  between  one-minimum  and  monotonic  ally  increasing  potentials — fe/ / 
assumes  its  maximum  (os  0.22),  whereas  pe  and  paw  assume  their  respective 
minima  (cs  0.59  and  us  0.98,  respectively).  Hence,  the  cold-electron-beam  ap¬ 
proximation  is  least  accurate,  though  still  acceptable  for  a  first  estimate,  for 
<5  m  0.405;  it  becomes  better  with  decreasing  or  increasing  4  because  in  both 
cases  the  velocity  spread  decreases  while  the  beam  velocity  increases. 

The  thin-sheath  approximation  made  throughout  this  work  requires  that 
the  axial  extensions  of  both  the  emitter  and  the  collector  sheaths  be  much 
smaller  than  that  of  the  plasma  region.  In  a  Q  machine,  this  condition  is 
usually  satisfied,  except  when  the  collector  bias  is  raised  to  excessively  high 
values.60,51  Quantitative  estimates  of  collector  sheath  widths  can  be  obtained, 
e.g.,  from  Fig.  3  of  Ref.  48. 


In  accordance  with  the  foregoing  considerations,  we  now  replace  a  long 
one-emitter  diode  of  the  above  kind  with  an  “approximate”  diode  having  in¬ 
finitely  thin  sheaths  and  a  uniform  plasma  region  of  length  L.  The  electrons 
of  the  approximate  diode  form  a  cold  beam  whose  plasma  density  and  velocity 


are  respectively  equal  to  the  plasma  density  and  average  velocity  of  the  orig¬ 
inal  diode.  Apart  from  the  non-trivial  external  circuits  we  have  in  mind,  this 
approximate  diode  is  formally  identical  with  the  uniform  classical  Pierce  diode 
described  in  Subeec.  II.A.  One  can  thus  associate  with  it  a  Pierce  parameter  a 
as  defined  by  Eq.  (1).  In  an  extension  of  the  concepts  introduced  in  Ref.  48,  we 
now  propose  a  method  yielding  quantitative  estimates  of  a.  This  establishes  a 
link  between  experimental  operating  conditions  and  theoretical  results  such  as 
given  in  Sec.  VI  and  elsewhere.33,34 


In  terms  of  standard  nondimensional  quantities, 43,40  Eq.  (1)  can  be  re¬ 
written  in  the  form 


where  Pp  =  flp/n,  is  the  normalised  plasma  density, 


(3) 


is  the  Richardson  emission  density,  A*  and  W*  are  the  effective  Richardson 
constant  and  work  function  for  thermal  electron  emission  from  the  emitter, 
A  —  L/l,  is  the  normalised  system  length,  and  l,  =  (xr/4irn(e3)1^3  is  the 
scaling  length.  With  the  numerical  results  displayed  in  Fig.  2b  of  Ref.  48  we 
have  calculated  (a/A)  as  a  function  of  the  neutralisation  parameter,  the  result 
being  shown  in  Fig.  2. 


For  a  given  experimental  situation,  the  Pierce  parameter  a  can  now  be 


obtained  as  follows.  With  the  method  proposed  in  Sec.  3  of  Ref.  48,  one  first 
determines  <5  and  A.  One  then  finds,  by  inspection  of  Fig.  2,  the  appropriate 
value  of  (<5/A),  which,  upon  multiplication  by  A,  yields  the  desired  result. 


In  order  to  get  a  feeling  for  the  numbers  that  may  occur  in  practice,  we 
find  it  useful  to  consider  yet  another  representation.  Multiplying  both  sides  of 
Eq.  (3)  by  I,,  we  obtain 


a 

I 


-£(*)■ 


(4) 


which  quantity  is  to  be  plotted  against  the  equilibrium  plasma  density.  Fig* 
ure  3  shows  several  curves  (a/L)  versus  ftp,  each  curve  corresponding  to  a 
specific  emitter  temperature.  The  emitter  material  assumed  is  tungsten  with 
a  Richardson  constant  A*  =  120  A  cm*3  K*3  and  effective  work  functions  for 
electron  emission  and  surface  ionisation  W*  =  W*  =  4.52  V,  cf.  Fig.  3a  of 
Ref.  47.  Each  curve  has  an  upper  and  a  lower  branch,  corresponding  to  one- 
minimum  and  monotonically  increasing  potentials,  respectively.  The  arrows 
indicate  the  direction  of  increasing  4.  These  curves  can  be  interpreted  as  fol¬ 
lows.  Consider  a  diode  with  L  =20  cm  and  T  =  2200  K,  and  assume  that 
the  method  of  Ref.  48  has  predicted  a  one-minimum  equilibrium  distribution 
with  ftp  » l.lxl07cm~s.  From  Fig.  3  we  find  (a/L  =  2.6cm-1,  and  hence 
a  =  52  —  16.6*.  Upon  solving  the  characteristic  equation  (Subsec.  V.A)  we 
would  find  that  there  are,  in  the  short-circuit  case,  two  unstable  eigenmodes 
with  frequencies  wi  =  (0.035  +  0.032* )wp  and  ui?  =  (0.13  +  0.0025t)wp.  However, 
anticipating  the  results  of  Subsec.  VI. B  we  expect  that  both  of  these  modes  can 
be  stabilised  by,  e.g.,  inserting  a  sufficiently  large  resistance  into  the  external 
circuit. 


Generally  speaking,  the  curves  plotted  in  Fig.  3  suggest  that,  for  typical 
operation  conditions,  the  Pierce  parameter  a  will  usually  lie  well  above  *,  so 
that,  in  the  short-circuit  case,  the  positively  biased  single-ended  Q  machine 
should  always  be  Pierce  unstable.  In  fact,  the  experimental  evidence  is  that 
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it  is  always  unstable41,42  except  for  very  high  collector  biases, 50,51  but  at  the 
present  stage  of  understanding  it  does  not  seem  possible  to  clearly  indicate  the 
dominant  instability  mechanism.  E.g.,  it  has  been  pointed  out4*  52  that  the  pos¬ 
itively  biased  single-ended  Q  machine  may  be  Buneman  unstable  as  well,  but 
the  results  there  were  obtained  in  the  infinite-plasma  approximation  and  there¬ 
fore  have  only  preliminary  character.  For  large  values  of  a,  the  growth  rates  of 
the  classical  Pierce  instability  become  much  smaller  than  the  plasma  frequency, 
cf.  Subsec.  VT.A.  They  can,  moreover,  be  modified  by  a  number  of  compo¬ 
nents  such  as  external-circuit  elements,  plasmarwall  interaction  processes,  finite 
sheath  widths,  collisions,  finite  velocity  spreads,  or  two-dimensional  geometry. 
If,  however,  a  is  not  too  far  above  r,  the  Pierce  instability  is  likely  to  manifest 
itself  more  clearly,  and  the  theoretical  results  given  in  Subsec.  VI.  may  not  be 
all  too  unrealistic. 

In  any  case,  we  believe  that  the  present  understanding  of  bounded  plasma 


systems — even  of  relatively  simple  ones — is  far  from  complete,  and  that  much 
more  theoretical  and  experimental  work  will  be  needed  to  improve  this  situation 
to  a  significant  extent. 


m.  DERIVING  THE  INTEGRAL  EQUATIONS 


A.  Specifications  (level-one  quantities) 

In  this  subsection  we  formally  define  the  perturbational  problem  to  be 
solved.  The  level-one  qnantities  as  listed  in  Table  I/I  are  specified  in  accordance 
with  the  general  description  of  the  problem  given  in  Secs.  I  and  II. 

We  restrict  ourselves  to  considering  “electron*  phenomena,  i.e.,  phenomena 
sufficiently  fast  that  the  ions  are  practically  unable  to  react  and  hence  can  be 
considered  to  be  immobile.  Thus,  the  electrons  are  the  only  species  we  have  to 
deal  with  (n„  =  1),  and  the  species  superscript  <x  can  be  dropped. 

In  the  plasma  region,  the  equilibrium  velocity  distribution  function  of  the 
electrons  is  assumed  to  be 

/(«)=*  ft,  *(v-0),  (5) 

which  corresponds  to  a  cold  electron  beam  of  density  ft,,  travelling  to  the  right 
with  velocity  0. 

Let  the  initial  perturbation  of  the  electron  distribution  function  be  of  the 

form 

fi(x,  v)  =  fti  6{x  -  i)  S(v  -  0) .  (6) 

This  means  that  the  final  solution  will,  among  other  things,  contain  the  “spa¬ 
tial”  Green’s  function  of  the  problem,  from  which  the  solution  for  any  arbitrary 
initial  perturbation  profile  can  be  easily  constructed  by  suitable  integration  over 
f,  cf.  Subsec.  III.B. 

The  boundary  conditions  for  the  particles  are  specified  by  requiring  con¬ 
stant  (non-sero)  emission  at  the  left-hand  electrode  (fi(v  >  0,  t)  =  0),  and  con- 


stant  (sero)  emission  at  the  right-hand  electrode  (fr(v  <  0,t)  =  0).  In  terms  of 
the  functions  introduced  in  Eq.  (26)/I,  these  conditions  can  be  expressed  as 


f,g(v  >  0,  t)  =  0 

(7°) 

bf*'  (w  >  0, «'  <  0)  =  0 

(76) 

&(«<  o,*)  =  o 

(7c) 

K°\v<  O.w'  >  0)  =0. 

(7  d) 

The  exemplary  external  circuit  we  consider  consists  of  a  d.c.  voltage  source 
shaping  the  equilibrium,  a  small-signal  generator  SG,  a  resistance  A,  and  an 
inductance  t,  all  connected  in  series  as  shown  in  Fig.  1.  (Clearly,  A  and  L 
are  defined  per  unit  cross-section  area.)  For  this  system,  the  external-circuit 
condition  (33)/I  combined  with  the  potential-balance  equation  (35)/I  reads 
A V,(t)  -  AVp(t) 

A  (8) 

=  V„S(t  -  r)  +  Aj e(t)  +  I— j,(0 , 

where  SG  has  been  assumed  such  as  to  produce  a  delta-function  signal  at  time 
r,  with  r  >  0.  Thus,  the  solution  will  also  contain  the  “temporal”  Green’s 
function,  from  which  one  may  ultimately  obtain  the  solution  for  any  arbitrary 
external  signal  Ve(r ),  cf.  Subsec.  m.B. 

In  order  to  completely  specify  the  initial  state  of  the  external  circuit  (and, 
hence,  of  the  whole  diode  system),  one  has  to  prescribe  either  two  of  the  three 
initial  values  ( dje/dt)i ,  and  AVj»,  whence  the  third  one  is  readily  obtained 

via  (8).  In  what  follows,  and  (dje/dt)i  will  be  assumed  to  be  given. 

Following  the  method  developped  in  Part  I,  the  initial-value  and  external- 
perturbation  problem  just  defined  must  be  solved  simultaneously  for  the  func¬ 
tions  J«(t),  E(x,t),  fi(v  >  0, t),  and  fr(v  <  0, t).  However,  since  we  have 
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prescribed  the  two  latter  functions  to  vanish  identically,  the  problem  actually 
reduces  to  finding  je  and  E. 


B.  Basic  modes  and  general  solution 

The  linear  perturbational  problem  specified  in  Subsec.  III.A  involves  the 
four  amplitude  factors  fa,  V„,  fa,  and  {dje/dt)i,  which  can  be  prescribed 
independently  of  each  other.  Hence  we  conclude  that  the  solution  will  be  a 
linear  superposition  of  four  independent  ‘basic*  modes,  which  we  choose  to 
define  as  shown  in  Table  I.  These  bask  modes  must  not  be  confused  with  the 
eigenmodes  discussed  in  Sec.  V;  they  rather  involve  sperpositions  of  eigenmodes, 
cf.  Sec.  IV.  Bask  modes  1  and  2  are  the  spatial  and  temporal  Green’s  functions, 
respectively,  whereas  bask  modes  3  and  4  are  connected  with  the  initial  state 
of  the  external  circuit. 

As  opposed  to  the  b ask-mode  problem  specified  in  Subsec.  HI.  A  and  solved 
in  Sec.  IV,  the  general  perturbational  problem  for  the  Pierce-type  configuration 
considered  is  obtained  by  replacing  (6)  with 

fi(*,  v)  =  fi(x)6(v  -  0)  (9) 


and  (8)  with 


A  ?.(*)  =  AVp(t) 

«?.(*)  +  + 


(10) 


where  the  initial  perturbation  profile  /,(x)  and  the  external-generator  signal 
Ve(t)  are  now  arbitrary.  The  solution  to  this  problem  can  be  constructed  by 
linear  superposition  of  basic  modes  in  the  form 


Thus,  solving  the  full  problem  is  practically  equivalent  to  finding  the  basic 
modes;  only  these  will  be  dealt  with  in  the  rest  of  the  present  section  and  in 
Sec.  IV. 


C.  Higher-level  quantities 


Starting  from  the  level-one  quantities  specified  in  Subsec.  III. A,  we  now 
construct  the  level-two  through  level-four  quantities  according  to  Table  I/I.  To 
this  end  it  proves  convenient  to  first  introduce  the  auxiliary  functions 

=  (- +  - >)  <«*> 
K,H  -  ,+!*  Y  4  fap .  , .  j2+fstt')  (12c) 

'  '  2t L  (w  +  0wp)3  V  •  0  /  '  ; 

n,  /(|ir»  WP°  XT'  1  Lr*(W+W(i“  0 


-1  -  » 


ci/  +  put* 


■(L-Z)) 


—  .  .  x  .WpO  \ — > 

itfilx,  w)  = - 7  +1—  >  7 - 77 

1  ’  L  2L  ££  (u  +  0wpy 


(ilw  +  0w„)x  ,  .w  +  /9wp  \ 

exp  — - — - 1  -  t - -  x  ) 

®  ®  / 


which  are  non-singular  for  all  values  of  w. 


(25) 


*<’•">  *  /.  {«•  [(- !)  *£=»  -  } 


+  Ve1 


ARi(*,w)1  f^Hi(*,w)l 

iL  u  j  ^  dt  J  |  xL  w  | 

(26) 

$t(v  >  0,  [x'],u/)  =  0 

(27) 

kio(v  >  0,w)  =  0 

(28) 

Sf(t»  <  0,(x'],w)  =  0 

(29) 

kis(v  <  0,u>)  =  0 

(30) 

(Remark:  The  operator*  Voi,  Vo,,  Vj,  and  Vr  are  not  required  here  because  the 
function*  upon  which  they  act  vanish  identically.) 


D.  The  integral  and  coefficient  equations 

With  the  level-four  quantities  calculated  in  Subsec.  m.C,  we  can  now  ex¬ 
plicitly  establish  the  integral  equations  (37)/I  through  (40)/I,  as  well  as  the 
equivalent  matrix  equation  (48)/I,  for  the  extended  Pierce  problem  defined  in 
Subsec.  m.A. 

Equation  (37)/I  remains  unchanged: 

+  fI‘dx,E(x,,»)  =  -V-rfM-  (31) 

Jo 

However,  due  to  the  simple  boundary  conditions  (7),  Eq.  (38)/l  reduces  to 

-k6(x,w)je(w)  +  E(x,io)  =  fcs(z,w).  (32) 


while  Eqs.  (39)/I  and  (40)/I  are  identically  satisfied  and  hence  need  not  be 
considered  in  what  follows. 


T 


We  thus  have  to  solve  Eqs.  (31)  and  (32)  for  je(w)  and  E(x,u).  The 
equivalent  matrix  equation  for  the  expansion  coefficients  is  given  by  Eq.  (48) /I: 


D(u)  u(o>)  =  k(w) ,  (33) 

where  the  coefficient  matrix  D(u),  the  known-coefficient  vector  k(w),  and  the 
unknown-coefficient  vector  u(u>)  are  now  considerably  reduced  in  comparison 
with  their  general  counterparts  (49a,b,c) /I: 


--(t) 


(34) 

i  / 

II 

Hr5 

(35a,  b) 

The  determinant  of  (34)  is  easily  found  in  the  analytical  form 

l£MI  «  ZM  + 


=  Ze(u)  +  f  dxk3(x,u) 

Jo 

=  (A  —  tut)  +  R3(u)  . 

(Jj 


IV.  SOLVING  FOR  THE  BASIC  MODES 


A.  Solution  of  the  integral  equations. 

The  eolation*  j,(u)  and  E(x,u>)  to  Eqa.  (31)  and  (32)  can  be  found  ana¬ 
lytically.  On  inserting  E(x,m),  as  calculated  from  (32),  into  (31),  and  taking 
into  account  (36),  we  readily  obtain 


AM 


-V’.oM  -  Jq  dxkt(x,ui) 

LCMI 


(37) 


Both  the  numerator  and  the  denominator  of  this  expression  contain  terms  pro¬ 
portional  to  w-1  and  thus  become  singular  for  u  — ►  0.  Since  these  singularities 
are  inessential  but  inconvenient  with  respect  to  what  follows,  we  multiply  both 
the  numerator  and  the  denominator  of  (37)  by  «.  With  a  view  to  the  basic 
modes  defined  in  Subsec.  HLB,  we  then  re-write  (37)  in  the  form 


where 


G(w)  s«|£(w)| 

=  w  (A-iw  I)  +  4in L  Rs(u) , 


(39) 


H(u>)  =  w  |  -V^w)  -  £  dx  ka(x, c)  J 


=  /«  ^(l,(w|C)  +  vtr  £,2,(w|r)  +  Jet  A(3'M 


(40) 
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^l*(wl€)  =  — 4irte 
A(2,Hr)  =  — w«*"T 

£|3|(u>)  =s  fa;£  +  t'A  A3(u>) 

£,4,(w)  =  »£  Jf3(w), 


u> 


(I-  -  £)  /?3 (<*^)  +  ^  —  ^("If) 
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jT’MO* 

li2](u\r)  m 

23,M  = 


G(w) 

£<2>(u>|r) 


G(u) 

£<3>M 


G(w) 


*  1  ’  G(w) 


Inserting  (38)  into  (32)  yields 


£(*,«)  =  +  *•(*»") 


=  /«^U,(*|W|0  +  ^r^,a,(*.wk) 


+  ^<3>(x,w)+(^j  ^4'(X,U,) 


where 


^^(x.wjC)  =  ^*(*IW)  +  *s1,(*.«IO 

#2,(x,«/|r)  =;J2,(w|r)fcs(x,u;) 

^,3,(x,w)  =  ^3,(w)fc8(x,w)  +  &g31  (x,  w) 
£,4,(x,w)  =  j<4,(u)k6(x,u)  +  ^g4'  (x,  w)  • 


(41a) 

(41b) 

(41c) 

(41d) 


(42a) 

(42b) 

(42c) 

(43d) 


(43) 


(44a) 

(44b) 

(44c) 

(44d) 


The  fnnctions  fcg1*  through  itg4*  can  be  found  in  a  straightforward  manner  by 
inspection  of  Gq.  (26). 
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With  (37)  and  (43)  we  have  found  the  solutions  to  Eqs.  (31)  and  (32),  the 
time  Laplace  transforms  of  the  basic  modes  being  given  by  Eqs.  (42)  and  (44). 

B.  Inverse  Laplace  transformation  and  time-dependent  basic  modes 

The  Laplace  transforms  jJm,(w)  and  ^m)(x,w)  (m  =  1,...,4),  as  given 
by  Eqs.  (42)  and  (44),  are  explicit  analytical  functions  of  w  and  can  be  easily 
inverted  into  the  time-dependent  basic  modes  jjm,(t)  and  £|m,(x,t),  cf.  Table 
I.  Since  in  the  present  example  the  Laplace  transforms  exhibit  no  branch  cuts, 
the  inversion  procedure  requires  the  evaluation  of  pole  contributions  only,  cf. 
Sec.  n.G/I  and  Fig.  2/1.  The  poles  to  be  considered  are  the  plasma  poles  (w  = 
±wp)  and  the  eigenfrequencies  (w  =  w„,  v  =  . . . ,  1, 2, . . .).  According  to  Eqs. 
(60) /I  and  (39),  the  latter  are  the  seros  of  the  determinant  (36)  or,  equivalently, 
of  the  function  G(w): 

G(w„ ;  L,  0,  ftp ,  A,  L)  =  0 ,  (45) 

where  all  parameters  are  shown  explicitly  for  later  purposes. 

After  some  algebra,  application  of  the  inverse  time  Laplace  transformation 
(8b) /I  to  (42)  and  (44)  leads  to  the  time-dependent  basic  modes 
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where  G1  =  cLG/du.  Since  the  asymptotic  part  in  (47b)  may  not  contain  a  con* 
at  ant  term,  we  conclude  that  =  0.  From  the  basic  modes  (46) 

and  (47),  the  general  solution  to  our  pert ur bat ional  problem  can  be  constructed 
as  described  in  Subeec.  ffl.B.  For  given  initial  perturbations  and  external- 
generator  signals,  this  is  a  purely  mechanical  task  which  will  not  be  pursued 
here  any  further.  Let  us  recall  that  jj1*  and  are  the  spatial  Green’s  func¬ 
tions  (for  initial  plasma  perturbations),  and  ji2>  and  are  the  temporal 
Green’s  functions  (for  external-generator  signals). 

By  inspection  of  (46)  we  see  that  the  external-circuit  current  densities 
Jim|(t)  (m  =  1, ...  ,4)  are  just  superpositions  of  exponentials  exp(-*o>„t),  i.e., 
their  time  behavior  is  exclusively  governed  by  the  eigenfrequencies  While 
bask  modes  1,  3,  and  4  are  due  to  initial  perturbations  and  hence  give  non-sero 
contributions  right  from  time  t  —  0  on,  basic  mode  2  is  activated  by  an  external 
delta  pulse  at  time  r  and  hence  contributes  only  for  t  >  r. 

According  to  (47),  the  corresponding  electric-field  patterns  are  much  more 
complex  in  that  each  of  them  exhibits  an  initial  transient  and  an  asymptotic 
solution.  At  a  given  position  z,  the  initial  transient  is  on  during  the  intervals 
0  <  t  <  z/0  (for  basic  modes  1,3,4)  and  r  <  t  <  r  +  x/0  (for  basic  mode  2). 
Its  time  behavior  is  not  only  governed  by  the  eigenfrequencies  but  also  by  the 
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plasma  frequency.  The  left-hand  boundary  of  the  transient  pattern  moves  to 
the  right  as  x  =  et  and  x  =  0(t  —  r),  respectively,  thus  giving  way  to  the  “dean” 
eigenmode  pattern  of  the  asymptotic  solution.  After  t  =  L/0  and  t  =  r  +  L/0, 
respectively,  the  transient  is  no  longer  present,  and  only  eigenmodes  are  left  in 
the  system.  Inspection  of  (47a-d)  shows  that  the  electric-field  profile  associated 
with  an  eigenfrequency  is  essentially  given  by  the  function  #1(2,  w„)  as 
defined  by  (12a). 

After  a  sufficiently  long  time,  the  ‘‘uppermost*  eigenmode — i.e.,  the  one 
with  the  most  positive  Imu^ — will  be  dominant. 
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V.  CHARACTERISTIC  EQUATION,  EIGENFREQUENCIES, 
AND  EIGENMODES 


With  the  basic-mode  solutions  (46)  and  (47),  the  detailed  linear  evolution 
of  our  Pierce-type  diode  can  be  calculated  for  arbitrary  initial  and  external  per¬ 
turbations  (9),  (10)  as  indicated  by  Eq.  (11).  The  full  details  of  this  evolution 
may  be  required,  e.g.,  for  interpreting  an  experiment  where  the  propagation 
of  externally  generated  small-amplitude  pulses  is  studied.  In  many  applica¬ 
tions,  however,  it  will  be  sufficient  to  concentrate  on  the  asymptotic  behavior. 
Whichever  is  the  case,  the  evaluation  of  Eqs.  (46)  and  (47)  requires  quantitative 
determination  of  the  eigenfrequencies  w„.  These  and  the  associated  eigenmode 
profiles  are  the  subject  of  the  present  section. 


A.  Normalised  characteristic  equation 

Quite  discour agingly,  the  eigenfrequencies  w„  are  functions  of  no  less  than 
five  physical  parameters,  namely  of  the  quantities  L,  0,  flp,  A,  and  t  appearing 
in  the  characteristic  equation,  or  ‘dispersion  relation”,  (45).  However,  the  com¬ 
plexity  of  the  problem  can  be  substantially  reduced  by  introducing  appropriate 
normalisations.  Let  us  re-write  Eq.  (45)  in  the  nondimensional  form 

d(riv-,a,p,\)  =  0  (48) 

where 

?  =  ^  (49) 

Up 

is  the  normalised  frequency,  a  is  the  classical  Pierce  parameter  defined  by  Eq. 

(1). 

(50  a) 


Rw2 

p=  ^ 

H  4x€ 
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are  the  normalised  external-circnit  resistance  and  inductance,  respectively,  and 


d'(ij;a,p,A) 


ojp  G(u) 

4xt0 


=  [rf  -  l)-2  {e,an  [(»73  +  l)  sin  a  +  2t*jcosa] 
— \tf*  —  ipi) 5  +  (a  +  2A)»74  +  2iprf  —  (a  +  A)t>2 


*(2  +  p)ri}  ■ 
(51) 


We  thus  have  to  solve  for  the  roots  r\u  of  Eq.  (48),  which  are  the  normal¬ 
ised  eigenfrequencies  and  “only*  depend  on  the  three  parameters  a,  p,  and  A. 
Each  of  these  parameters  may  vary  between  0  and  oo,  and  for  each  parame¬ 
ter  triple  there  are  infinitely  many  eigenfrequencies  rj„,  whose  imaginary  parts 
are  bounded  from  above.  The  question  thus  arises  of  how  to  investigate  in  a 
systematic  manner  the  relevant  eigenfrequencies  and  the  associated  eigenmode 
profiles  in  a  given  parameter  domain  of  interest. 


B.  Eigenfrequencies 

For  a  given  parameter  triple  (a,  p,A),  all  eigenfrequencies  down  to  some 
chosen  (Imt})m«n  can  be  approximately  located  as  follows.  One  starts  out 
by  plotting  the  complex  image  curves  (j(ry,  a,  p,  A)  for  a  suitable  sequence  of 
complex  straight-line  portions  of  the  form  (-M  <  Ret)  <  M,  Im»j  =  const.), 
where  M  is  positive  and  sufficiently  large  to  simulate  infinity,  and  1m  ij  assumes 
successively  lower  values  from  some  suitable  (Imr;)m<u  down  to  (Im  r;)min.  By 
inspection  of  the  image  curves  one  has  to  find  out  which  ones  come  close  to  the 
origin  G  =  0,  and  at  which  values  of  Re  rj.  One  thus  obtains  the  starting  values 
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needed  for  the  numerical  eolation  of  Eq.  (48).  This  procedure  is  closely  related 
to  the  classical  Nyquist  method53  and  can,  in  principle,  be  automated. 


Due  to  the  symmetry 

<5(-Rer;  +  tImn;a,p,A)  =  [<5(+Req +t  Imn;a,P,A)]  ,  (52) 

which  is  not  difficult  to  verify  on  Eq.  (51),  we  have  that  eigenfrequencies  with 
a  non-vanishing  real  part  always  occur  in  pairs  of  the  form 

^^iRe^'+tlm  (53) 

where  the  superscript  (d)  indicates  the  “doublet"  character  of  these  eigenfre¬ 
quencies,  and  Re  oj*4*  is  defined  to  be  positive.  Similarly,  we  will  label  purely 
imaginary,  or  "singlet",  roots  of  Eq.  (48)  by  the  superscript  (a). 

As  will  become  apparent  in  the  following  subsection,  both  eigenfrequencies 
of  a  doublet  (53)  correspond  to  one  and  the  same  eigenmode,  so  that  it  is  conve¬ 
nient  to  distinguish  between  the  "eigenfrequency  index"  v  and  the  "eigenmode 
index”  p.  Thus,  a  doublet  (53)  corresponds  to  two  values  of  v  but  to  just  one 
value  of  p. 

C.  Egenmode  profiles 

By  inspecting  the  asymptotic  parts  of  the  basic-mode  solutions  (47a-d)  we 
see  that  the  complex  electric  field  associated  with  an  eigenfrequency  =  upr\u 
is  generally  given  by 

4c'(*,t)  =CJJ,(i,ai,)riM+l1  (54) 

with  C  and  7  real.  This  equation  can  be  re-written  in  the  form 
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Elei(x,t)  =  C{Ref?i(z,w„)  coa(Rewl/t  +  -7) 

+  ImAifz,^)  8in(Rewl/t  +  7)}  e1™"-* 

f  IT  (55) 

+  iC  | Re  Ri(x,u/V)  cos(ReWi,t  +  7  +  — ) 

+  sin(Reul/t  +  7  +  ^)}  elmu',t , 

from  which  it  is  evident  that  the  real  and  imaginary  parts  describe  the  same 
spatial  and  temporal  behavior,  except  for  an  amplitude  factor  *  and  a  phase 
difference  x/2.  Hence,  we  may,  without  loss  of  generality,  define  a  real  and 
normalised  eigenmode  profile  by  the  real  part  of  Eq.  (55),  with  C  —  1  and 

7=0: 

JE £"*’  =  {Re  f?i (z,u„)  cos(Rew„t) 

(w) 

+  Im J2i(z,(i>t,)  sin(Rew„t)} 

We  furthermore  choose  to  define  the  associated  potential  profile  in  the  normal¬ 
ised  form 


vyn]  (x,  =  <**'  Krn)  (*'.  *) 


=  {ReR6(x,u„)  coe(Rew„t)  +  ImA((z,uv)  sin(Rewut)}  elmUl' 
Since  the  functions  Ri(x,u)  and  Rs{x,ui)  exhibit  the  symmetry 
Ri.s(x,  ~  Rew  +  ilmu)  =  [Ri,5(zi  +  Rew  + 1  Imw)f  , 


Eqs.  (56)  and  (57)  imply 

e{;;]  (z,  t)  =  (z,  t)  =  E';n)(x,  t)  (59a) 

vj;nl  (*,  t)  =  V™  (z,  t)  =  Vf"  (z,  t) ,  (596) 

so  that  both  eigenfrequencies  of  a  doublet  (53)  correspond  to  one  and  the  same 
physical  eigenmode. 
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Clearly,  Re  Ri(x,wv)  and  Re  ^(x,  w„)  represent  the  initial  values  of  the 
eigenmode  profiles  (56)  and  (57),  respectively.  One  easily  verifies  that 
ReRi(0,uv)  =  1  and  ImRi(0,uv)  =  ReRs(0, w„)  =  ImRefO.ui,,)  —  0.  Since 
the  constants  mission  assumption  (7a)  implies  zero  space-charge  perturbation 
at  x  =  xt  —  0,  we  also  expect  3(Ref?i)/5z  =  3(Im Ry)/dx  =  0  for  x  — ►  0. 
For  a  singlet  eigenfrequency  (Rewj,*1  =  0),  both  Im  Ri  and  Im  Rf,  vanish  iden¬ 
tically,  and  the  eigenmodes  are  either  purely  damped  (ImwJ,*'  <  0)  or  purely 
growing  (Imo/j,*1  >  0),  thus  retaining  their  initial  spatial  profiles.  (We  have 
not  investigated  the  details  of  the  singular  case  =  0.)  For  doublet  eigenfre- 
quencies  (Rewj,d)  /  0),  Im  J2i  (z,ujtd*)  and  Im/26(x,u/ j,**’)  represent  the  spatial 
field  and  potential  structures  after  a  quarter-period,  R ewj,d)t  =  x/2,  and  the 
eigenmode  is  a  superposition  of  two  damped  (Imwji^  <  0),  constant-amplitude 
(Imwj»d|  —  0),  or  growing  (Imwj,4’  >  0)  standing  waves  with  different  spa¬ 
tial  profiles  and  a  temporal  phase  difference  of  |,  respectively  resulting  in  a 
damped,  constant-amplitude,  or  growing  propagating-wave-like  structure. 

In  view  of  the  foregoing  remarks,  the  spatial  field  and  potential  structures 
of  an  eigenmode  corresponding  to  a  singlet  eigenfrequency  are  uniquely 
represented  by  Re  •Ri(x,wj,*1)  and  Re  f?6(x,  wj/1),  respectively.  Representation 
of  a  doublet  eigenmode,  on  the  other  hand,  in  addition  requires  the  functions 
ImAi(x,c<;jid>)  and  ImR&(x,u>l,d>). 


VL  NUMERICAL  RESULTS 


In  this  section  we  present  systematic  numerical  results  for  the  uppermost 
normalised  eigenfrequencies  and  the  associated  eigenmode  profiles  in  some  typi-  I 

cal  parameter  domains,  cf.  Sec.  V.  In  Subsec.  VI. A.,  the  classical,  short-circuited 
Pierce  diode,  which  lends  itself  as  a  convenient  standard  case,  is  reviewed, 

a  classification  scheme  for  the  eigenfrequencies  is  introduced,  and  some  new  ^ 

eigenfrequencies  are  calculated.  The  effects  of  external-circuit  resistance  and 
external-circuit  inductance  on  the  eigenfrequencies  and  eigenmode  profiles  are 
considered  in  Subsecs.  VI. B  and  VI.C,  respectively.  Some  general  aspects  and 

implications  are  discussed  in  Subsec.  VI. D.  ^ 

A.  The  classical,  short-circuited  Pierce  diode  (a  >  o,  p  =  o,  A  =  o) 

The  classical  Pierce  diode,3  defined  in  Subsec.  II.A,  offers  itself  as  a  refer-  I 

ence  case  to  which  all  modified  configurations  can  be  conveniently  compared. 

Numerical  results  for  the  dominant  eigenfrequencies  have  been  given  before, 33,34 
but  for  the  systematic  description  attempted  here  we  have  found  it  necessary  to 
complement  these,  in  the  parameter  domains  considered,  with  new  results  for 
some  of  the  lower  (with  respect  to  Im  ij )  eigenfrequencies.  Moreover,  since  only 
little  information  on  eigenmode  structures  can  be  found  in  the  literature,13  17 
we  present,  throughout  this  section,  some  typical  eigenmode  profiles.  Follow¬ 
ing  Subsec.  V.C,  these  are  represented  in  terms  of  the  functions  ReAi,  Im  Ri, 

ReRs,  and  ImiZg. 


Throughout  the  present  section,  we  restrict  ourselves  to  varying  the  Pierce 
parameter  in  the  range  0.16jt  <  a  <  3 ir,  for  which  the  uppermost  eigenfre¬ 
quencies  are  shown  in  Fig.  4.  This  range  is  representative  in  that  it  covers 


the  “archetypical”  a-intervals  0  <  a  <  r,  r  <  a  <  2x,  and  2r  <  a  <  3», 
which  will  be  seen  to  exhibit  all  possible  modes  of  linear  behavior.  Since  these 
three  intervals  are  distinctly  different  from  each  other,  it  is  convenient  to  num¬ 
ber  the  eigenfrequencies  uniformly  within  each  of  them.  For  a  given  a-interval 
nr  <  a  <  (n  +  l)w,  where  n  =  0, 1, 2, . . .,  we  choose  to  number  the  eigenfre- 
quencies  %  so  as  to  reflect  the  ordering  of  Im  in  the  middle  of  the  interval, 
Le.,  at  o  =  (b  +  j)r.  Thus,  the  eigenfrequency  singlet  or  doublet  with  the 
uppermost  imaginary  part  there  is  assigned  p  =  1  throughout  the  interval,  the 
singlet  or  doublet  with  the  next  lower  imaginary  part  there  is  assigned  ft  —  2 
throughout  the  interval,  etc.  If  /t  =  2  corresponds  to  say  an  eigenfrequency 
doublet,  we  will  speak  of  “eigenmode  2(<<l  ”,  etc. 

Throughout  the  three  a-intervals  shown  in  Fig.  4,  the  ordering  of  the 
eigenfrequencies  2(<<l  and  3<4>  remains  the  same,  so  that  no  re-numbering  is 
necessary  at  the  transition  between  intervals.  Above  3<**> ,  there  exist  no  other 
eigenfrequencies  than  the  ones  shown.  With  increasing  a,  the  packing  den¬ 
sity  of  the  eigenfrequencies  with  respect  to  Imq  increases  rapidly.  Thus,  at 
a  =  0.5*  =  1.5708  there  are  only  five  eigenfrequencies  with  Im»j  >  -3.5  (a 
doublet  counts  for  two  eigenfrequencies),  at  a  =  1.5*  —  4.7124  there  are  six 
eigenfrequencies  with  Imr;  >  -1,  whereas  at  a  =  2.5r  =  7.8540  we  already 
find  44  eigenfrequencies  with  Imrj  >  —  1.  Eigenfrequencies  lying  below  3^)  are 
not  shown,  with  the  exception  of  the  singlets  41**  (in  r  <  a  <  2 r)  and  X*'1  (in 
2r  <  a  <  3x;  order  ft  unknown). 

Figure  4  shows  that  the  dominant  eigenmodes  (n  =  1),  which  determine  the 
linear  time-asymptotic  behavior,  are  non-oscillatory  stable,  or  purely  damped, 
for  0  <  a  <  r,  non-oscillatory  unstable,  or  purely  growing,  for  r  <  a  <  2r, 
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oscillatory  unstable  for  2x  <  a  <  3t  —  e  (where  e  ~  0.3),  and  both  oscillatory 
and  non-oscillatory  stable  for  3jt  —  e  <  a  <  3 x.  The  pattern  of  the  last  two 
intervals  basically  repeats  with  increasing  a,  so  that  for  3x  <  a  <  4ir  we  have 
non-oscillatory  instability,  etc.;  for  large  values  of  a,  the  growth  rates  decrease 
approximately  as  a-1.33,34 

Electric-field  and  electric-potential  profiles  of  the  dominant  eigenmodea  for 
a  »  O.Sx,  1.5 r,  and  2.5*  are  shown  in  Figs.  6  and  12,  Figs.  8  and  14,  and 
Figs.  10  and  16,  respectively.  For  the  classical  case,  the  field  and  potential 
structures  become  more  complex  with  increasing  a,  and  the  potential  profiles 
always  return  to  sero  at  x  —  L. 

B.  The  purely  resistive  case  (a  >  o,p  >  0,  X  =  o) 

The  purpose  of  the  present  subsection  is  to  establish  how  the  eigenfre- 
quencies  and  eigenmode  structures  change  if  we  allow  for  finite  resistance,  but 
still  require  sero  inductance,  in  the  external  circuit.  We  decide  to  answer  this 
question  for  the  three  a-values  1.5v,  1.5ir,  and  2.5 w,  as  representatives  of  the 
a-intervals  they  are  located  in,  cf.  Subsec.  VI. A.  For  each  a,  the  imaginary  and 
real  parts  of  the  normalised  eigenfrequencies  rj„  are  plotted  for  normalised  re¬ 
sistances  in  the  range  0  <  p  <  100,  and  some  typical  field  and  potential  profiles 
are  given  for  the  dominant  eigenmodes.  The  value  p  =  0  retrieves  the  classical 
short-circuit  case,  while  p  — ►  oo  corresponds  to  the  open-circuit,  or  floating, 
diode. 

Fig.  5  shows  the  eigenfrequencies  r^*1,  and  for  a  =  0.5 ir.  No  other 
eigenfrequencies  were  found  above  3(<i>.  We  observe  that  increasing  resistance 
results  in  weaker  damping  for  l*'1  but  stronger  damping  for  2,<,,  and  3,d|.  Thus, 


the  system  becomes  “less  stable*  with  increasing  resistance,  and  tends  towards 
marginal  stability  for  large  p.  Figure  6  shows  that  the  eigenmode  profiles  of 
l<d*  change  significantly  with  p. 

The  eigenfrequencies  »;['*,  rjj*11,  and  for  a  =  1.5x  are  shown  in 
Fig.  7.  No  other  eigenfreqnencies  were  found  above  4**1  and  3(d),  respectively. 
The  most  remarkable  result  is  that  l**1  remains  unstable.  However,  Imi)[** 
decreases  monotonkally  with  increasing  p,  and  the  system  approaches  marginal 
stability  as  p  — *  oo.  Figure  7  also  shows  that  ImfjJ*1  crosses  Imt?^*  at  p  ~  3, 
so  that  from  there  on  4**1  is  the  third  highest  mode.  Profiles  of  eigenmode  l1** 
for  different  values  of  p  are  given  in  Fig.  8. 

Figure  9  shows  the  eigenfreqnencies  of  the  modes  l|d),  2<d|,  3,d),  X'**, 
X*'1',  and  X(,>"  for  a  =  2.5 r.  No  other  eigenfreqnencies  were  found  above 
3,d|.  Here  we  encounter  a  much  more  complex  and  interesting  behavior  than 
for  a  ~  0.5tr  (Fig.  5)  and  a  —  1.5*  (Fig.  7).  The  classically  unstable  mode 
l(d>  is  stabilised  at  p  as  4.  At  p  ~  19.2,  ri[d)  is  "absorbed*  by  a  singlet  mode, 
which  is  seen  to  be  the  continuation  of  the  classically  high-order  mode  X1*1, 
cf.  Fig.  4.  Quite  strangely,  the  imaginary  part  of  describes  an  S-shaped 
curve  which,  for  19.2  <  p  <  93,  exhibits  three  branches  denoted  by  X(*’, 
and  X'*1".  At  p  ~  93,  a  doublet  branches  off  to  the  right,  which  we  interpret 
as  the  continuation  of  l<d|.  Returning  to  the  uppermost  eigenfrequencies  we 
may  state  that  the  dominating  eigenmode  is  oscillatory  growing  for  0  <  p  <  4, 
oscillatory  damped  for  4  <  p  <  19.2,  and  purely  damped  for  p  >  19.2.  For 
large  p,  the  system  tends  towards  marginal  instability.  Eigenmode  profiles  of 
l(dl  and  X’*'"  for  different  values  of  p  are  given  in  Fig.  10. 


C.  The  purely  inductive  case  (a  >  o,  p  =  0,  A  >  o) 

In  this  subsection  we  investigate  the  eigenfrequencies  and  eigenmode  struc¬ 
tures  for  the  case  of  finite  inductance  but  tero  resistance  in  the  external  circuit. 
As  in  Subsec.  VLB,  we  restrict  ourselves  to  the  three  a-values  0.5*,  1.5*,  and 
2.5*.  For  each  a,  the  uppermost  eigenfrequencies  are  plotted  in  the  range 
0  <  A  <  100,  and  some  typical  eigenmode  profiles  are  given.  The  value  A  =  0 
recovers  the  classical  short-circuit  case  of  Subsec.  VI. A,  whereas  A  — *  oo  corre¬ 
sponds  to  a  current-stabilised  external  circuit. 

Fig.  11  shows  the  dependence  on  A  of  the  classical  eigenfrequencies 
qJ,d|,  and  for  a  =  0.5*.  However,  we  also  observe  a  new  doublet  L|d| 
(eigenfrequency  rj^),  which  does  not  exist  for  A  =  0  and  thus  is  inherently 
due  to  the  finite  inductance  in  the  external  circuit.  Since  Imq^d>  always  lies 
well  above  the  imaginary  parts  of  the  other  eigenfrequencies  shown — no  others 
were  found  above  3(d|  — the  “inductive*  mode  L,d|  is  the  dominant  one  and 
even  goes  unstabV  within  some  small  intervals  of  the  region  0  <  A  <  0.09.  This 
means  that  the  finite  inductance  now  present  has  changed  the  linear  response 
behavior  fundamentally  as  compared  with  the  short-circuit  case,  which  is  non- 
oscillatory  stable.  Fig.  12  shows  some  eigenmode  profiles  of  Lfd*  for  different 
values  of  A.  For  A  — ►  0,  the  oscillation  frequency  Re  tends  towards  infinity — 

which  behavior  is  also  observed  for  the  other  a-values  considered — and  thus 
precludes  numerical  solution  of  Eq.  (48)  for  very  small  values  of  A.  As  A  — ►  oo, 
r j £d)  approaches  sero,  whereas  the  damping  rates  of  the  other  modes  increase 
logarithmically.  With  Eqs.  (46a-d)  in  mind,  this  is  the  qualitative  behavior  we 
expect  for  the  case  of  a  current-stabilised  external  circuit. 


The  variation  with  A  of  the  eigenfrequencies  r)^d\  ijj1*,  r^'1,  and 
for  a  =  1.5t  is  shown  in  Fig.  13.  No  other  modes  were  found  above  3(d>. 
Contrary  to  a  =  0.5x,  the  classically  dominant  mode  l*'1  remains  dominant, 
retaining  its  purely  growing  character.  Its  field  and  potential  profiles  change 
with  A  as  shown  in  Fig.  14.  As  A  — *  oo,  both  1  and  r^* 1  tend  to  sero,  whereas 
the  damping  rates  of  the  other  eigenfrequencies  increase  logarithmically.  In 
some  intervals  of  the  region  A  <  1.1,  the  inductive  mode  L(d)  goes  unstable. 

Figure  15  shows  the  variation  with  A  of  the  eigenfrequencies  rj[d\  , 
t) 3d’,  and  rj^  for  a  =  2.5*.  No  other  eigenfrequencies  were  found  above  3,d*. 
The  stability  behavior  is  very  similar  to  that  for  a  =  1.5s-  (Fig.  13),  except  that 
now  the  dominant  mode  is  an  oscillatory  one.  Figure  16  shows  the  eigenmode 
profiles  of  l*d|  for  various  values  of  A.  The  inductive  mode  L|d|  introduces  ad¬ 
ditional  instability  within  various  intervals  in  the  region  A  <  4,  at  an  oscillation 
frequency  much  higher  than  that  of  l!d). 

D.  Discussion 

The  results  of  the  present  section  have  shown  that  finite  resistance  or  in¬ 
ductance  in  the  external  circuit  can  strongly  modify  the  eigenfrequencies  and 
eigenmode  profiles  as  compared  with  their  classical  counterparts. 

In  Subsec.  VI.B  we  have  found  that  increasing  the  external-circuit  resis¬ 
tance  can  change  the  dominant  mode  (i)  from  non-oscillatory  stable  to  non- 
oscillatory  marginally  stable  (a  =  0.5x),  (ii)  from  non-oscillatory  unstable  to 
non-oscillatory  marginally  unstable  (a  =  1.5x),  and  (iii)  from  oscillatory  unsta¬ 
ble  successively  to  oscillatory  stable,  non-oscillatory  stable,  and  non-oscillatory 
marginally  stable  (a  =  2.5*). 


Perhaps  the  most  interesting  finding  of  this  section  is  that  finite  external- 
circuit  inductance  gives  rise  to  a  new  oscillatory  mode  L,d*,  which  does  not 
exist  for  A  =  0.  The  occurrence  of  this  inductive  mode  is  plausible  because 
the  system  consisting  of  the  diode,  the  resistance  A,  and  the  inductance  L 
(Fig.  1)  may  be  looked  upon  as  an  RLC  circuit,  with  the  diode  acting  as  a 
capacitor  whose  properties  are  partly  determined  by  the  plasma  in  it.  Thus, 
the  eigenmode  L1***  may  be  interpreted  as  a  modified  RLC  oscillation.  For 
all  values  of  a  considered  here,  we  have  made  the  remarkable  observation  that 
there  are  regions  of  A  where  L,d’  is  unstable,  and  that  with  increasing  a  these 
regions  become  larger.  For  a  =  0.5 r,  L,d*  is  the  dominant  mode,  whereas  for 
a  =  1.5 ir  and  a  =  2.5*  it  is  second,  at  least  in  some  range  of  A. 

Globally  speaking,  we  have  found  that  the  linear  response  behavior  of  our 
Pierce-type  diode  depends  crucially  on  all  three  parameters  considered  here, 
namely  a,  p,  and  A.  This  induces  us  to  speculate  that  also  the  nonlinear  final 
states  may  strongly  depend  on  these  parameters,  as  well  as  on  others  that  may 
have  to  be  defined  for  more  realistic  cnfigurations.  These  final  states  can  be 
time-independent,  oscillatory,  or  turbulent,34  and  it  appears  that  much  more 
work  is  needed  to  clarify  these  phenomena,  the  parameter  ranges  in  which 
they  occur,  and  their  dependencies  on  the  parameters  themselves.  Particle 
simulations  tackling  some  of  these  problems  are  under  way  and  will  be  reported 
elsewhere.43 

It  should  be  borne  in  mind  that  the  three  cases  considered  in  Subsecs.  VI. A, 
VI.B,  and  VI.C  are  highly  idealised  ones,  thus  singling  out  specific  external- 
circuit  effects.  Realistic  external  circuits,  on  the  other  hand,  always  involve 
resistive,  inductive,  and  capacitive  components  in  some  more  or  less  complex 


Apart  from  deriving  these  detailed  results,  we  have  also  attempted  to  give 
a  general  account  of  the  Pierce  problem  and  its  implications.  In  so  doing  it 
has  been  observed  that  the  Pierce  problem  is  potentially  relevant  to  a  wide 
class  of  laboratory  plasmas  and  has  recently  been  attracting  a  remarkable  deal 
of  interest  in  the  context  of  various  modern  concepts  and  applications.  More 
specifically,  we  have  argued  that  the  low-density  single-ended  Q  machine  can, 
to  a  certain  extent,  be  approximated  by  a  simple  Pierce-type  configuration. 

In  spite  of  the  fairly  large  number  of  papers  that  have  already  been  written 
on  Pierce-type  problems,  we  have  to  conclude  that  only  rather  restricted  aspects 
have  been  studied  so  far,  and  that  substantial  progress  in  the  theory  of  bounded 
plasma  systems  will  require  more  comprehensive  and  realistic  approaches.  A 
first  step  in  this  direction  has  been  proposed  in  Part  I,  but  major  extensions — 
e.g.,  to  more  general  geometries — are  still  needed. 
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TABLE  I.  Definition  of  basic  modes. 
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FIG.  1.  Schematic  of  extended  Pierce-type  configuration. 

FIG.  2.  Pierce  parameter  over  normalised  system  length  versus  neutralisation 
parameter. 

FIG.  3.  Pierce  parameter  over  system  length  versus  equilibrium  plasma  den¬ 
sity. 

FIG.  4.  Imaginary  parts  (a)  and  real  parts  (b)  of  the  normalised  eigenfrequen- 
cies  i 7^’,  ijj4*,  t and  for  [0.16*r  <  a  <  3x|p  =  0|A  =  0]. 

FIG.  5.  Imaginary  parts  (a)  and  real  parts  (b)  of  the  normalised  eigenfrequen- 
cies  r^d\  and  for  (a  =  0.5t  |0  <  p  <  100 |  A  =  0]. 

FIG.  6.  Electrostatic  field  (a)  and  potential  (b)  profiles  of  eigenmode  l*'*  for 
[o«0.5jr|p  =  0,l,  10, 100 1 A  ~  Oj. 

FIG.  7.  Imaginary  parts  (a)  and  real  parts  (b)  of  the  normalised  eigenfrequen- 
ciea  i jj*',  and  for  [a  *  1.5jt  1 0  <  p  <  100 1 A  =  0]. 

FIG.  8.  Electrostatic  field  (a)  and  potential  (b)  profiles  of  eigenmode  ll#l  for 
[a«  1.5*|p  =  0, 1,10, 100 1 A  =0]. 

FIG.  9.  Imaginary  parts  (a)  and  real  parts  (b)  of  the  normalised  eigen&equen- 
ciea  t)[d\  ijg*',  X***',  and  for  [a  =  2.5x|0  <  p  <  100 1  A  =  0). 

FIG.  10.  Electrostatic  field  (a,b)  and  potential  (c,d)  profiles  of  eigenmode  l(d| 
for  [a  =  2.5*  |  p  =  0, 4, 19 1 A  =  0],  and  of  eigenmode  for  [a  =  2.5ir  |  p  = 
30, 100 1 A  =  0]. 

FIG.  11.  Imaginary  parts  (a)  and  real  parts  (b)  of  the  normalised  eigenfre- 
quencies  r?1/1,  rf^d\  ifg**,  and  for  [a  =  0.5ir  |  p  =  0 1 0  <  A  <  100]. 
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FIG.  12.  Electrostatic  field  (a,b)  and  potential  (c,d)  profiles  of  eigenmode  Lldi 
for  [a  =  0.5*  |  p  =  0 1 A  =  0.035, 0.07, 1, 10, 100]. 

FIG.  13.  Imaginary  parts  (a)  and  real  parts  (b)  of  the  normalised  eigenfre- 
qnencies  ni*1.  ^ ,  *»i‘\  and  ^  for  fa  =  1.5ir|p  =  0 }  0  <  A  <  100]. 

FIG.  14.  Electroetatk  field  (a)  and  potential  (b)  profiles  of  eigenmode  l1*1  for 
[a  =  1.5 1  p  =  0 1 A  =  0, 1,10,100]. 

FIG.  15.  Imaginary  parts  (a)  and  real  parts  (b)  of  the  eigenfreqencies 
rtf',  i £°,  and  for  [a  =  2.5*  |  p  =  0 1 0  <  A  <  100]. 

FIG.  16.  Electrostatic  field  (a,b)  and  potential  (c,d)  profiles  of  eigenmode  1(<1> 
for  [a  *  2.5*  |  p  =  0 1 A  =  0, 1, 10, 100]. 


